Measurement of direct photon emission in 
K + — > 7r + 7r°7 decay using stopped positive 

kaons 



M.A. Aliev a , Y. Asano b , T. Baker , P. Depommier d , 
M. HasinofP, K. Horie ,:L , Y. Igarashi c , J. Imazato , 
A.P. Ivashkin a , M.M. Khabibullin a , A.N. Khotjantsev a , 
Y.G. Kudenko a , A.S. Levchenko a , G.Y. Lim c , 
J. A. Macdonald 8 ', O.V. Mineev a , C. Rangacharyulu h , 
S. Sawada c and S. Shimizu f '* 

KEK-E470 Collaboration 

a Institute for Nuclear Research, Russian Academy of Sciences, Moscow 117312, 

Russia 

b Institute of Applied Physics, University of Tsukuba, Ibaraki 305-0006, Japan 

c Institute of Particle and Nuclear Studies (IPNS), High Energy Accelerator 
Research Organization (KEK), Ibaraki 305-0801, Japan 

d Laboratoire de Physique Nucleaire, Universite de Montreal, Montreal, Quebec, 

Canada H3C 3J7 

e Department of Physics and Astronomy, University of British Columbia, 
Vancouver, Canada V6T 1Z1 

f Department of Physics, Osaka University, Osaka 560-0043, Japan 

g TRIUMF, Vancouver, British Columbia, Canada V6T 2A3 

h Department of Physics, University of Saskatchewan, Saskatoon, Canada S7N 

5E2 



Abstract 

The radiative decay K + — ► tt + -k j (K^^) has been measured with stopped positive 
kaons. A K W 2~ sample containing 4k events was analyzed, and the K n 2y branching 
ratio of the direct photon emission process was determined to be [6.1 ± 2.5(stat) ± 
1.9(syst)] x 10 -6 . No interference pattern with internal bremsstrahlung was observed. 
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1 Introduction 



The chiral anomaly [1], which is a basic feature of quantum field theory, has 
been the subject of extensive theoretical investigations. There are a number of 
calculations for decays of the pseudoscalar mesons - n, 77, and K [2]; however 
few experimental studies of this anomaly have been performed so far. With the 
increased interest in chiral perturbation theory (ChPT) [3], the chiral anomaly 
has attracted renewed attention. In particular, non-leptonic kaon decays are 
well-suited for the study of ChPT. It is known, however, that the genuine 
manifestation of the chiral anomaly in non-leptonic decays is restricted to the 
radiative decay of K + — > 7r + 7r°7 (K w2l ) for the charged kaon [4]. 

The K + — > 7r + 7r {K %2 ) decay is hindered because it violates the AI—1/2 
rule. As a consequence, the internal bremsstrahlung (IB) contribution to K n2l 
decay is suppressed, although it is still dominant. This feature, in turn, could 
enhance the direct photon emission (DE) which is sensitive to meson structure. 
Here, the chiral anomaly appears as the magnetic component of DE. In terms 
of ChPT, it enters the magnetic amplitude at 0(p i ), while IB arises from 
0(p 2 ) [4]. Also, the electric amplitude of DE may arise from 0(p 4 ), which 
one can identify through the interference pattern (INT) with IB. Thus the 
K n2l decay channel is one of the most interesting and important channels in 
determining the low energy structure of QCD. 

Three in-flight-if 4 " experiments [5,6,7] and one stopped-.fr 4 " experiment [8] 
have reported on the DE branching ratio in the 7r + kinetic energy region of 
55 < T n + < 90 MeV. These experiments restricted the T n + region to avoid 
backgrounds from K + — > 7r + 7r°7r° (K n3 ) decay, although the DE component 
is significantly enhanced with decreasing n + energy. The weighted average of 
the in-flight K + experiments is Br(K w2l , DE) = (18 ± 4) x 1(T 6 [9] which is 
much larger than the result from the stopped K + experiment Br(K n2j , DE) = 
(4.7±0.8±0.3) x 10 -6 [8]. Therefore, further experimental studies are necessary 
to settle this discrepancy of the DE branching ratio. It is also desirable to check 
the existence of the INT component which has not yet been observed. 

In this letter, we present a new measurement of direct photon emission in 
K n 2y decay. The experiment used a stopped K + beam in conjunction with 
a superconducting toroidal spectrometer with a 12-sector iron-core [10]. We 
measured the K^ 2l events by extending the T n + region (T w + >35 MeV) below 
the K„ 3 upper threshold (T^ ax =55 MeV). 
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Fig. 1. Cross sectional side view of the E470 setup. Assembly detail of the active 
target, TOF1, and ring counters is shown in the inset. 

2 Experiment 



The experiment was performed at the KEK 12 GeV proton synchrotron. The 
experimental apparatus was based on the E246 experiment searching for T- 
violating transverse muon polarization in K + — > 7i°n + v (K^) decay [11]. Be- 
sides the T- violation search, spectroscopic studies for K + — > 7r°e + z/ (K e3 ), 
K^, and K n3 decays have been carried out using the same detector sys- 
tem [12,13,14]. A schematic cross sectional view of the detector is shown in 
Fig. 1. Because of the rotational symmetry of the 12 identical gaps in the 
spectrometer and the large directional acceptance of the photon detector [15], 
distortions of the K n2 -y spectra due to instrumental misalignment, if any, were 
drastically reduced. 

A separated 660 MeV/c K + beam (ti + /K + ~ 7) was used. Each K + was 
discriminated from pion background by a Cerenkov counter with an efficiency 
of more than 99%. The typical K + beam intensity was 1.5 x 10 5 over a 2 s 
spill duration with a 4 s repetition rate. The kaons were slowed down by a 
degrader and stopped in an active target located at the center of the detec- 
tor system. The K n2 f events were identified by analyzing the n + momentum 
with the spectrometer and detecting three photons in the CsI(Tl) calorimeter. 
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Charged particles from the target were tracked and momentum-analyzed us- 
ing multi-wire proportional chambers at the entrance (C2) and exit (C3 and 
C4) of the magnet gap, as well as an array of ring counters [16] surround- 
ing the active target system. The n + s were discriminated from e + and /i + by 
calculating the mass squared (M| OF ) of the charged particles from the time- 
of-flight between TOF1 and TOF2 counters. Counter TOF1 surrounded the 
active target and counter TOF2 was located at the exit of the spectrome- 
ter. The measurement was carried out for the central magnetic field strength 
B=0.65 T which was chosen to match the spectrometer acceptance to the ir + 
momentum distribution due to the DE process. 

The photon detector, an assembly of 768 CsI(Tl) crystals, covered 75% of the 
total solid angle. There were 12 holes for outgoing charged particles to enter 
the spectrometer and 2 holes for beam entrance and exit, as shown in Fig. 1. 
The photon energy and hit position were obtained, respectively, by summing 
the energy deposits and taking the energy-weighted centroid of the crystals 
sharing a shower. Timing information from each module was used to identify a 
photon cluster and to suppress accidental backgrounds due to beam particles. 
A Pb-plastic sandwich detector with 2.6 radiation length was set at the outer 
radius of the magnet pole to monitor photons passing through the photon 
detector holes, as shown in Fig. 1. 



3 K n2 -y event selection 

K n 2-y decays at rest were obtained by the following procedure which was care- 
fully studied in a Monte Carlo simulation (see below). The K + decay time, 
defined as the n + signal at the TOF1 counter, was required to be more than 
1.4 ns later than the K + arrival time measured by the Cerenkov counter. 
This reduced the fraction of K + decays in-flight contamination to the level 
of 10 -3 . Events with n + decays in-flight and scattering of charged particles 
from the magnet pole faces were eliminated by requiring the hit position in 
ring counters to be consistent with the charged particle track. The selection 
of 7T+, by 15000 < M| OF < 30000 MeV 2 /c 4 , rejected e+s and /i+s. To remove 
K W 2 events, the reconstructed n + momentum p n + (Fig. 2(a)), corrected for 
the energy loss in the target, was required to be p n + <175 MeV/c. Events 
with three clusters in the CsI(Tl) calorimeter were selected by assigning two 
photons from 7r° — > 71 + 72 and a radiated photon (73), while events with other 
cluster number were rejected. The opening angle between the two photons was 
required to be larger than 37° to eliminate photon conversion backgrounds in 
the active target system. Also, the photon conversion events could be removed 
by selecting only events in which there was one charged particle hit in TOF1 
as well as in ring counters. The total K + energy M K + = E v + + Y%=i anc ^ 
momentum vector p K + = p n + + Y%=i P 7i were used to reject background con- 
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taminations. The windows of 420 < M K + < 500 MeV/c 2 and \(p K +)x, y ,z\ <50 
MeV/c were imposed. 

The final step in reconstructing K^ 2l was to determine which photons were a 
pair from a 7r°, because there are three possible combinations to form the 7r°. 
The IB events with incorrect pairing mimic the DE decay and might introduce 
a systematic error. In order to find the correct pairing, a quantity Q 2 defined 

as, 



Q 2 = (M 7172 - M n0 y/<j 2 M o + (E w0 + E rA + E n+ - M k+ ) 2 /<j 2 Mk+ 
+ (cos9^-cos9^-5 1 ) 2 /al +y3 
+ (cose A . L o C " cos« S - S 2 ) 2 /al +7r0 , 



was calculated for the three combinations, and the pair with the minimum 
Q 2 , <5min> was adopted as the n° pair. Here, M 7l72 is the invariant mass of the 
selected pair and E^o is the n° energy obtained as, 



E n0 = ^E\ + 2£ 7l £ 72 cos# 7l72 + E\ + m 2 , (2) 

where # 7l72 is the opening angle of the selected pair and m n o is the rest mass of 
7T°. The angles 9 n + 73 and 9 w + n o are that between n + and 73 and that between 
7T + and 7T°, respectively. The superscripts MEAS and CALC stand for the 
measured angles and the calculated angles using the K n2l kinematics as, 



qCALC 



< , XLr m 2 K+ + m 2 + - m 2 - 2m K +E 7T + - 2m K +E J3 + 2E 13 E^ 
cos ^+73 =— 9 i„ I,- 1 1 1 > ( 3 ) 



/iCAlc m K+ + m l+ + m l" - 2rn K +E n + - 2m K +E 7T o + 2E 7T oE 7V + 
cosV n+7r0 = — , (4) 

A \Pn+\\Pn°\ 

where m, E, and p denote the rest mass, energy, and momentum vector, 
respectively, of the particle shown in the subscript. The a and offset val- 
ues for each term were om — H-4 MeV/c 2 , o"M if+ =10.5 MeV/c 2 , <r 7r + 73 =0.48, 
0-^0=0.06, M^. = 116.0 MeV/c 2 , M K = 471.6 MeV/c 2 , ^=0.35, and 5 2 =0.06. 
The choice of Eq. (1) and these parameters was determined to get the highest 
probability for the correct pairing by using a simulation with the IB compo- 
nent. The correct pairing probability was estimated to be 85% for IB, 79% 
for DE, and 83% for INT. Since K w3 and K n2 decays do not satisfy the K n2l 
kinematics, the Qmin values of the K n3 and K w2 events are essentially larger, 
as shown in Fig. 2(b). Therefore, the cut of <5 2 nin < 4.3 drastically reduced 
these contaminations. 



The number of good K n2 ^ events is 4434 after the above standard selection 
conditions. The background fraction of the K nS decay was estimated to be 
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P jt+ (MeV/c) P (MeV/c) 



Fig. 2. (a) p n + distribution and (b) correlation plot of p n + and Q^in without p n + 
and Q^in cuts. To demonstrate background rejection performance by the Q m i n cut, 
the p w + cut is not applied here. 

1.2% from the simulation. The solid lines in Fig. 3 show the experimental K n2l 
spectra after subtracting the backgrounds: (a) 7r + momentum, (b) opening 
angle between n + and 73, (c) opening angle between tt° and 73, and (d) 73 
energy. These are compared with the calculated ones (see below). 



4 Monte Carlo simulation 



In order to obtain the detector response function, a Monte Carlo simulation 
was carried out for both the 7r + measurement by the spectrometer and the 
photon measurement by the CsI(Tl) calorimeter using a GEANT-based Monte 
Carlo code. In this simulation, K n2l events due to the IB, DE, and INT pro- 
cesses were generated with the corresponding matrix elements [17] and an- 
alyzed in the same manner as the experimental data. The accuracy of the 
simulation was carefully checked by comparing the K n2 spectrum with the 
experimental one. The decay was also used to check the reproducibility 
of the experimental conditions in the simulation. 

The detector acceptance ratios were obtained to be f2(DE)/fi(IB) = 2.63 and 
f2(INT)/f2(IB) = 2.56 which were used to determine relative decay strength of 
the DE and INT components normalized to the IB component. Here it should 
be noted that events with incorrect pairing were included in the sample, 
and this effect was also treated in the simulation. The pairing probability 
could be controlled by requiring extra cuts in addition to the standard cuts at 
the cost of good K n2l events. Two cut conditions, (I) Q 2 CC — Q^m > 1 an d (H) 
cos# 7172 > or cos0„-o 73 < 0, were individually or simultaneously imposed, 
where Q 2 CC is the second smallest Q 2 value among the three combinations. 
The pairing probability and good event acceptance under these conditions are 
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Fig. 3. K^-y spectra: (a) ir + momentum, (b) opening angle between n + and 73, 
(c) opening angle between 7r° and 73, (d) 73 energy. The solid lines are the experi- 
mental data. The dotted and dashed lines show the fitted spectra and the sole DE 
components, respectively, under the assumption of no INT component. The bump 
structure in (b) is due to the CsI(Tl) assembly structure with 12 holes. 

Table 1 

Summary of the detector acceptance O and correct pairing probability P for IB and 
DE. 



Condition 


ft(IB) 


ft(DE) 


P(IB) 


P(DE) 




(ID" 3 ) 


(10- 3 ) 


(%) 


(%) 


standard 


0.24 


0.63 


85 


79 


standard® (I) 


0.20 


0.53 


91 


84 


standard® (II) 


0.22 


0.53 


92 


80 


standard® (I)® (II) 


0.18 


0.46 


95 


85 



summarized in Table 1. It is possible to increase the probability up to 95% 
(IB) and 85% (DE), although 23% of the K n2 f events are rejected. These cuts 
were used for the estimation of systematic errors. 
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Fig. 4. Q^ in distributions of (a) and (b) K n3 events for the experimental data 
(solid line) and the simulation (dotted line). Dashed line in (a) shows the simulated 
K-xz background survived after the K n 2y selection. 

The background was also estimated by the simulation and subtracted 
from the experimental data. Since the present spectrometer acceptance over- 
lapped the K n3 region, the estimation of the K n3 fraction was important. The 
obvious K n3 events were identified as 3 photon clusters with one escaping 
photon using cuts of T n + < 55 MeV (p w + < 135 MeV/c), 350 < M K + < 440 
MeV/c 2 , 15000 < M| OF < 35000 MeV 2 /c 4 , and Q 2 lin > 10. These events 
were compared with the simulation to confirm the accuracy of the simulation 
and to check the correctness of the background estimation. Fig. 4(a), (b) show 
the <5min distributions of the and K n3 events, respectively, for the experi- 
mental data (solid line) and the simulation (dotted line). They are in excellent 
agreement, supporting the correct estimation of the K w3 fraction in the 
sample. The dashed line in Fig. 4(a) shows the simulated K n3 background 
which survived after the K n2 ^ selection. 



5 Results 

The DE branching ratio was determined using the three observables, cos6 , 7r + 73 , 
cos0„-o 73 , and E ri which have characteristic spectra corresponding to the de- 
cay processes. The experimental distribution p(cos8 n + ri , cos9 n o 13 , E l3 ) of 
data from the standard cut conditions was fitted to a simulation spectrum of 
IB[l+a(DE/IB)+/3(INT/IB)] with a and (3 being free parameters. The DE 
and INT components were obtained to be a=(0.92±°;|j)% and /3=(-0.58la!a)%, 
respectively. The reduced x 2 is 1-27. Fig. 5 shows the x 2 contour plot in 
the (a, (3) space. (3 is consistent with zero and we conclude that there is 
no significant INT component observed within the experimental accuracy. 
By assuming (3 = in the fitting, the DE component was determined to 
be a=(0.76±0.31)% with a reduced x 2 =l-28. The DE branching ratio was 



8 



5.D 




-5.0 



_L 



_L 



1.0 2.0 
<x(%) 



3.0 



Fig. 5. x 2 contour plot in the (a, f3) space. Contour levels of 1-, 2-, and 3-a 
constraints are drawn. 
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Fig. 6. (a) cos6' 7r + 73 and (b) cos^o^ spectra of the DE component obtained by 
subtracting the fit IB component from the experimental data (dots) and the fitted 
DE component (histogram). 

derived to be BriK^^ DE) = (6.1±2.5) xl0~ 6 by normalizing it to the the- 
oretical value of the IB branching ratio [17]. The dotted and dashed lines in 
Fig. 3 show the fitted spectra and the sole DE components, respectively, under 
the assumption of (3 =0. It can be seen that 26% of DE comes from the region 
T n + < 55 MeV. Fig. 6(a), (b) show the cos8 n + 13 and cos^o^ spectra, respec- 
tively, of the DE component obtained by subtracting the fitted IB component 
from the experimental data (dots) and the fitted DE component (histogram). 

The major systematic errors in the DE measurement come mainly from imper- 
fect reproducibility of the experimental conditions in the simulation as follows. 
The 7 mispairing effect was studied by varying the pairing probability. In ad- 
dition to the standard K n2l selection conditions, extra cuts given in Table. 1 
were imposed. The DE branching ratio was found to be consistent within 
Aa <0.17% which was regarded as the upper bound of this systematic error. 
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Table 2 

Summary of estimated systematic uncertainty. All items are added in quadrature 
to get the total. 

Error source Uncertainty of a (%) 

7 mispairing <0.17 

7 calibration 0.13 

K W 3 background 0.03 

Accidental background 0.08 

Spectrometer field ambiguity 0.02 

Total systematic error <0.23 

The uncertainty due to the calibration error of the CsI(Tl) calorimeter was 
evaluated by introducing the maximum conceivable errors in the simulation, 
and found to be Aa = 0.13%. The K w3 fraction in the K^ 2l data depends 
slightly on the K n3 decay form factors, however this effect is negligibly small. 
The systematic error due to the ambiguity of the fraction was obtained 
as Aa = 0.03% by rejecting events with T n <55 MeV. The contribution from 
the accidental background was estimated to be Aa = 0.08% by changing its 
fraction varying the TDC gate widths. The uncertainty due to the spectrom- 
eter field ambiguity was obtained to be Aa = 0.02%. As shown in Table 2, 
the total systematic error was estimated to be Aa to tai <0.23% by adding all 
the item in quadrature. However, this upper bound was regarded here as one 
standard deviation of the systematic error. 



6 Conclusion 

We have performed a new measurement of the K + — > 7r + 7r°7 decay by ex- 
tending the T w + lower bound down to 35 MeV. The detector response and 
acceptance functions were evaluated by a Monte Carlo simulation based on 
GEANT. The reproducibility of the experimental conditions in the simula- 
tion was carefully checked using K n2 and K n ^ decays. The background 
contamination in the K n2 -y data sample was also estimated by the simulation. 

Comparing the experimental data with the simulation, the DE component was 
obtained to be 

a = [0.76 ± 0.31(stat) ± 0.23(syst)]%, 

under the assumption of no INT component. The DE branching ratio was then 
determined to be 
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Br(DE, total) = [6.1 ± 2.5(stat) ± 1.9(syst)] x 10~ 6 , 



by normalizing it to the theoretical value of the IB branching ratio [17]. In 
order to compare this result with the previous ones, the partial branching ratio 
in the region of 55 < T w + < 90 MeV was also calculated as 



Br(DE, partial) = [3.2 ± 1.3(stat) ± l.O(syst)] x 1(T 6 , 

which is consistent with the previous stopped experiment [8]. We did not 
observe any significant effect due to the INT component, indicating the pure 
magnetic nature of the direct photon emission process in decay. 
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